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We study the power of upcoming weak lensing surveys to probe dark energy. Dark energy mod- 
ifies the distance-redshift relation as well as the matter power spectrum, both of which affect 
the weak lensing convergence power spectrum. Some dark-energy models predict additional 
clustering on very large scales, but this probably cannot be detected by weak lensing alone due 
to cosmic variance. With reasonable prior information on other cosmological parameters, we 
find that a survey covering f 000 sq. deg. down to a limiting magnitude of R = 27 can impose 
constraints comparable to those expected from upcoming type la supernova and number-count 
surveys. This result, however, is contingent on the control of both observational and theoretical 
systematics. Concentrating on the latter, we find that the nonlinear power spectrum of matter 
perturbations and the redshift distribution of source galaxies both need to be determined accu- 
rately in order for weak lensing to achieve its full potential. Finally, we discuss the sensitivity 
of the three-point statistics to dark energy. 



I. INTRODUCTION 

Recent direct evidence for acceleration of the universe 
jl], ^] has spurred considerable activity in finding ways 
to probe the source of this acceleration, dark energy 
H [§ H §| (f° r a review of dark energy see Ref. 0]). Be- 
cause dark energy varies with redshift more slowly than 
matter, it starts contributing significantly to the expan- 
sion of the universe only relatively recently, at z < 2. 
This component is believed to be smooth (or nearly so), 
and therefore detectable mainly through its effect on the 
expansion rate of the universe. For these reasons, it 
is generally believed that type la supernovae (SNe la) 
and number-count surveys of galaxies and galaxy clus- 
ters have the most leverage to probe dark energy, as 
they probe the distance and volume in the desired red- 
shift range ||. Indeed, planned supernova surveys (e.g., 
SNAP[f] ) and number-count methods || are ex- 
pected to impose tight constraints on the smooth com- 
ponent; for example, <r(w) « 0.05 from SNAP, assuming 
a flat universe. 

The program of weak gravitational lensing (WL) is pri- 
marily oriented toward mapping the distribution of mat- 
ter in the universe. The paths of photons emitted by 
distant objects and travelling toward us are perturbed 
due to the intervening mass. The weak lensing regime 
corresponds to the intervening surface density of matter 
being much smaller than some critical value; in that case 
the observed objects (e.g. galaxies) are slightly distorted. 
The weak lensing distortions are small (roughly at the 1% 



[f ]http: / /snap. lbl.gov 



level) and one needs a large sample of foreground galaxies 
in order to separate the lensing effect from the "noise" 
represented by random orientations of galaxies. There- 
fore, observations of lensed galaxies provides information 
on the matter distribution in the universe, as well as the 
growth of density perturbations. Although the potential 
of WL has been recognized for around two decades (e.g. 
p"T[), only in the 1990s was there a surge of interest in 
this area |T^, |l4|, [f5|, 16, |f7j A unique property of 
WL is that it is sensitive directly to the amount of mass 
in the universe, avoiding the thorny issue of galaxy-to- 
mass bias. By measuring ellipticities of a large number 
of galaxies, one can in principle directly reconstruct the 
mass density field of an intervening massive object |fl8| . 
Indeed, the mass reconstruction of galaxy clusters has 
been successfully performed on a number of clusters (for 
a review, see Ref. Q). 

An exciting recent development, relevant to this work, 
was the discovery of weak lensing by large-scale struc- 
ture, announced by four groups [|o[ |2l], p3[ . The 
results are in mutual agreement and consistent with the- 
oretical expectations, which is remarkable given that they 
were obtained independently. Although current data im- 
pose weak constraints on cosmology (e.g., rule out the 
Einstein-deSitter Universe with f2jvf = I), future surveys 
with larger sky coverage and improved systematics are 
expected to impose interesting constraints on cosmolog- 
ical parameters. 

The goal of this work is to assess the power of weak 
lensing to constrain dark energy. This analysis therefore 
complements that of Huterer & Turner 0, where the 
efficacy of SNe la and number-count surveys was con- 
sidered. We follow the standard practice of considering 
dark energy to be a smooth component parameterized by 
its energy density (scaled to critical) fix and equation- 
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of-state ratio w = p/p [p3. Dark energy modifies the 
WL observables by altering the distance-redshift relation 
and the gr owth of density perturbations. As discussed in 
Sec. VI B, the nonlinear evolution of perturbations also 



depends on dark energy; this dependence is much more 
difficult to calculate and needs to be calibrated from N- 
body simulations. Overall, the dependence of WL on 
dark energy is somewhat indirect and expected to be 
weak, especially when degeneracy with other cosmolog- 
ical parameters is taken into account. Nevertheless, we 
shall show that, provided systematic errors are controlled 
and theoretical predictions sharpened, WL surveys can 
be efficient probes of dark energy, comparable to SNe 
la and number-counts. Proposed deep wide-field sur- 
veys such as LSST[f] , the aforementioned SNAP, and 
VISTA [|] will attempt to constrain dark energy through 
their WL programs, making our analysis particularly 
timely. 

Previous work on parameter determination from WL 
centered mostly on CIm and <7g, the rms density fluc- 
tuation in spheres of 8h _1 Mpc 16, 



here Hq 



100 h km/s/Mpc is the Hubble parameter today. Hu and 
Tegmark pq], however, used the Fisher matrix formal- 
ism to account for all 8 parameters upon which WL de- 
pends, and assumed dark energy to be the vacuum energy 
(therefore, fixed w = —1). We use the same set of param- 
eters, with two changes: we add w, and, guided by the 
ever-stronger evidence from the cosmic microwave back- 
ground (e.g., |27], |2^, |2^, |3^]), we assume a flat universe. 
To assess the accuracy of parameter determination, we 
too use the Fisher matrix machinery, which has proven 
to be an extremely efficient and accurate way to forecast 
errors in experiments where observables depend on many 
parameters. 

This paper is organized as follows. In Sec. |l| we go 
over the basic formalism and define the notation. In Sec- 



tions |Hj and |IV| we concentrate on the convergence power 
spectrum, and discuss its dependence on dark energy. 
Section ^ discusses the power of weak lensing surveys to 
probe dark energy, while Sec. [v| addresses systematic er- 
rors that can lead to biases in parameter estimation. In 
Sec. [VI I| we discuss the dependence of three-point statis- 
tics — bispectrum and skewness of the c onvergence — 
on dark energy. We conclude in Sec. VIII . 



II. PRELIMINARIES 

In this Section we cover the basic formalism of weak 
gravitational lensing (for detailed reviews, see Refs. |ll], 
pl[). We work in the Newtonian Gauge, where the per- 
turbed Friedmann-Robertson- Walker metric reads 



[fjhttp:/ /dmtelescope.org 
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ds 2 = - (1 + 2$) dt 2 + a 2 (t) (1 - 2$) x 
[d X 2 +r 2 {d6 2 + sin 2 0d(j) 2 )] 



(1) 



where we have set c — 1, x is the radial distance, $ is 
the gravitational potential, and k — 1, 0, —1 for closed, 
flat and open geometry respectively. We also use the 
coordinate distance r which is defined as 



r {x) 



{ (-K)-^ 2 s^h[{-Kf/ 2 X l if r» TOT < 1; 

x, if ^tot = l; 

i^- 1 /2 s i I1 (i^i/2 x ) ) if Otot > 1. 



where K is the curvature, SItot is the total energy den- 
sity relative to critical, and K = (Oxot — 1)^0- 

Gravitational lensing produces distortions of images 
of background galaxies. These distortions can be de- 
scribed as mapping between the source plane (S) and 
image plane (I) |32j 



Sxf = AtjSxj 



(3) 



where <5x are the displacement vectors in the two planes 
and A is the distortion matrix 
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FIG. 1: The assumed source galaxy distribution n(z). 



The deformation is described by the convergence k and 
complex shear (71,72). We are interested in the weak 
lensing limit, where |k|, \j\ <C 1. The convergence in any 
particular direction on the sky n is given by the integral 
along the line-of-sight 
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K(n, X )= / W( X ')S(x')d x ' 



(5) 



where S is the relative perturbation in matter energy den- 
sity and 



W( X ) = -n M H$g( X )(l + z) 
is referred to as the weight function. Furthermore 



9(X) = r(x) I d X 'n(x)^ }x V ' 



(6) 



r(x') 



r(x) r (Xs ~ X) 
r(Xs) 



(7) 

(8) 



where n(x) is the distribution of source galaxies in red- 
shift (normalized so that J dzn(z) = 1) and the second 
line holds only if all sources are at a single redshift z s . 
We use the distribution |Eol 



equation (see Eq. (||)), and is convenient for comparison 
with theory. However, it is very difficult to measure the 
convergence itself, as convergence depends on the mag- 
nification of galaxies which would somehow need to be 
measured (although there may be ways to do this; see 
Ref. J3^| ) [f ] . Note also that computing the convergence 
from the measured shear is difficult in general, since the 
inversion kernel is broad and requires knowledge of shear 
everywhere fL8[ . In the weak lensing limit, however, the 
problem is much easier, since the two-point correlation 
functions of shear and convergence are identical. 

In this work we use power spectrum of the convergence 
(defined in Eq. ( |l3|) below) as the principal observable 
that will convey information from weak lensing. 



III. CONVERGENCE POWER SPECTRUM 

The convergence can be transformed into multipole 
space (see e.g. Ref. 31 1) 



Kir, 



J dhn(h,x)Yj t 7 
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-z/zo 



(9) 



with zq = 0.5, which peaks at 2zq — 1 and is shown in 
Fig. [TJ. Our results depend very weakly on the shape of 
the distribution of source galaxies (assuming this distri- 
bution is known, of course). In particular, if all source 
galaxies are assumed to be at z = 1, the parameter un- 
certainties change by at most ~ 30% percent. Similarly, 
the distribution given by Eq. (|^) which peaks at z = 1.5 
would improve the parameter constraints by 20% or less. 

Some clarification is needed regarding observability vs. 
theoretical computability of WL quantities. The quantity 
that is most easily determined from observations is shear, 
which is directly related to the ellipticities of observed 
galaxies (in the weak lensing limit, shear is equal to the 
average ellipticity). Shear is given by Jl5| 



7i +ij2 = ~ {ip,xi 



- iip ■ 



(10) 



where ip is the projected Newtonian potential, ip,ij = 
—2 J g(x) dx, and commas denote derivatives with 
respect to directions perpendicular to the line of sight. 
Unfortunately, this quantity is not easily related to the 
distribution of matter in the universe and the cosmologi- 
cal parameters. Convergence, on the other hand, is given 
by 



The power spectrum of the convergence P ; K is then 
defined by 



Using Limber's approximation — the fact that the 
weight function W is much broader than the physical 
scale on which the perturbation 5 varies — the conver- 
gence power spectrum can be written as 



fx* w 2 (x) 



r 2 {x) 



(12) 



He 
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dz 



W 2 (z)r(z) A2 



H(z) 



A\l/r(z),z) (13) 



where in the second line we assume a flat universe where 
dx = dr. Here P(k, z) is the matter power spectrum as 
a function of redshift z, and 



A 2 (M) 



k 3 P(k,z) 
2^2 



(14) 



is power per unit logarithmic interval in wavenumber, 
which we also refer to as the matter power spectrum. 



\ ("0,11 +^,22) 



(11) 



which (in Limber's approximation) can be directly re- 
lated to the distribution of matter through the Poisson 



[f] There are two competing effects due to magnification of galax- 
ies: 1) "Magnification bias", the increase in the observed number 
of galaxies due to the fact that fainter ones can now be observed, 
and 2) increase in the apparent observed area on the sky due to 
lensing, which decreases the observed number density of galaxies. 
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FIG. 2: Top panel: The convergence power spectrum for 
three pairs of (fix, w). The shaded region represents 
l-cr uncertainties (corresponding to Qx = 0.7, w = — 1 
curve) plotted at each I. The uncertainties at low / are 
dominated by cosmic variance, and those at high I by 
Poisson (shot) noise; see Eq. (|l5|). We also show the con- 
tribution to Pi from the linear matter power spectrum 
only. Bottom Panel: Pf / APf ("signal-to-noise") for the 
convergence power spectrum for each individual I. 



1/2 

ered by a survey of dimension 9 and (7^) ~ 0.4 is 
the intrinsic ellipticity of galaxies. The first term cor- 
responds to cosmic variance which dominates on large 
scales, and the second to Poisson noise which arises due to 
small number of galaxies on small scales. Bottom panel 
of Fig. 2 shows the signal-to-noise Pf/APf. It is appar- 
ent that the bulk of cosmological constraints comes from 
multipoles between several hundred and several thou- 
sand. Wider and deeper surveys widen the range of scales 
with high signal-to-noise. Note also that the weak lens- 
ing power spectrum is relatively featureless because of 
the radial projection (Eq. (|l3|)). It can be characterized 
by amplitude (normalization), overall tilt, a "turnover" 
at I ~ 100 which is due to the turnover in the matter 
power spectrum, and a further increase at / ~ 1000 and 
flattening at / ~ 10000 which are due to the nonlinear 
clustering of matter. 

The systematic error in shear measurements ideally 
needs to be small enough so as not to exceed the sta- 
tistical error shown in Fig. ^. The maximum allowed 
systematic error can be estimated using the following ar- 
gument [f] . The shear variance in circular aperture of 
opening angle can be written in terms of the conver- 
gence power spectrum as fl3~i| 

(7 » (e)) = 2 .p^(AE) 2 {16) 




which is the power in unit logarithmic interval evaluated 
at I = 1/6 (here V l K = 1(1 + l)/(2ir)P l K ). The tightest 
requirement is on scales of 1 arcmin (I ss 2000), where 
the fractional uncertainty in power per unit logarithmic 
interval is about 1/200. Therefore the RMS of shear on 
scale 6 is given by 



Power spectrum of the convergence is displayed in the 
top panel of Fig. ^ for three values of (fix, w) and down 
to scales of about one arcminute (I = 10000). The un- 
certainty in the observed weak lensing spectrum is given 

by P|[t] 



AP? = 



(2£+l)f sky 



PP 



(15) 



where / s k y = 9 7r/129600 is the fraction of the sky cov- 



[f] Strictly speaking, Eq, ^ holds for Gaussian convergence field 
only. However, the non-Gaussianity of the convergence is milder 
than that of the matter due to the radia l projection which makes 
this a good approximation, see Sec. [VI ( 



1=1/6 



1=1/6 



= -MV. 

v ' is, 

4 XyiFIx 20() 
w 2.5 x 10 -5 

< 0.01 v^))- 



(19) 
(20) 

(21) 
(22) 
(23) 



Therefore, the systematic error in individual shear 
measurements should be less than 1% in order to be 
subdominant to statistical error — a very challenging 
requirement indeed! 



[fJProposed by M. Turner. 
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IV. DEPENDENCE ON DARK ENERGY 

The sensitivity of the convergence power spectrum to 
dark energy can be divided into two parts. Dark energy 

a) modifies the background evolution of the universe, 
and consequently the geometric factor W 2 {z)r(z) / H(z), 
and 

b) modifies the matter power spectrum. 
We now discuss each of these dependencies. 



A The lensing weight function 
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FIG. 3: The weight function W 2 (z)r(z) / 'H(z) for three 
pairs of (fix, w). 



Function W(z) is bell-shaped, and has a maximum at 
z w z s /2, where z s is redshift of lensed galaxies, indicat- 
ing that lensing is the most effective at distances halfway 
between the source and the observer. Since r(z) and H (z) 
are varying with redshift monotonically and slowly, the 
function W 2 (z)r(z) / H(z) will also be bell-shaped with 
maximum at z > z s /2. W(z), r(z) and 1/H(z) all de- 
crease with increasing w, and therefore the total weight 
decreases. As fix decreases, r(z) and l/H(z) decrease 
but W(z) increases, and the latter prevails; see Fig. [|. 
Therefore, changing w (fix) makes the normalization and 
total weight change with the same (opposite) sign, lead- 
ing to large (small) change in on large scales; see 
Fig. | 



B The Matter Power Spectrum 

The matter power spectrum can be written as 

A 2 (fc, z) = 8% + T 2 (fc, z) T NL (k, z\2A) 



where Sh is perturbation on Hubble scale today, T(k) is 
the transfer function, D(z)/D(0) is the growth of pertur- 
bations in linear theory relative to today, and T/vx(fc, z) 
is the prescription for nonlinear evolution of the power 
spectrum. 

In the presence of dark energy, the matter power spec- 
trum will be modified as follows. 

• The normalization Sh increases with increasing fix 
and decreasing w. This happens because the growth of 
structure is suppressed in the presence of dark energy, 
and the observed structure today can only be explained 
by a larger initial amount of perturbation. We choose to 
normalize the results to COBE measurements [34 , and 
adopt the fit to COBE data of Ma et al. (j35j, heretofore 
Ma QCDM) 

S H = 2 x 10- 5 a(0)- 1 fi^ +Caln(nAr) 

x exp [c 3 (n - 1) + c 4 (n - l) 2 ] (25) 



where C1...4 and cxq are functions of fix and w and are 
given in Ma QCDM. Since the COBE normalization for 
ACDM models is accurate to between 7% and 9% |3^, [37| , 
we adopt the accuracy of 10% for the dark-energy case. 

• The transfer function for cosmological models with 
neutrinos and the cosmological constant is given by fits 
of Hu and Eisenstein |38| , which we adopt in our analy- 
sis. These formulae are accurate to a few percent for the 
currently favored cosmology with low baryon abundance. 
Dark energy will not directly modify the transfer func- 
tion, except possibly on the largest observable scales (of 
size ~ Hq 1 ), where dark energy may cluster. This sig- 
nature can be ignored, as it shows up on scales too large 
to be probed by WL; we further discuss this in Sec. VG. 

• The linear theory growth function D(z) = S(z)/5(Q) 
can be computed from the fitting formula for the dark- 
energy models given in Ma QCDM, which generalizes the 
ACDM growth function formula of Carroll et al. 39 1 . We 
use this fitting function, noting that its high accuracy 
(< 2%) justifies avoiding the alternative of repeatedly 
evaluating the exact expression for the growth function 
(e.g. @, p. 341). 

• The last, and most uncertain, piece of the puzzle 
is the prescription for the non-linear evolution of den- 
sity perturbations. This is given by the recipe of Hamil- 
ton ct al. |^l[ as implemented by Peacock and Dodds 
(Q; heretofore PD), as well as Ma (Q, heretofore 
Ma ACDM). These two prescriptions were calibrated for 
ACDM models, although the PD formula seems to ade- 
quately fit models with w > — 1 (M. White, private com- 
munication). Ma QCDM prescription p5[ , on the other 
hand, gives explicit formulae for the nonlinear power 
spectrum in the presence of dark energy (i.e., a com- 
ponent with w > —1). Unfortunately, we found that the 
PD and Ma QCDM prescriptions agree (to ~ 15%) only 
at values of w where Ma QCDM was tested. At other 
values of w the maximum disagreement between the two 
is up to 50%, and it is not clear which fitting function, 
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if any, is to be used. We choose to use the PD prescrip- 
tion primarily because it is implemented for all w and 
therefore facilitates taking the derivativ e with respect to 
w needed for the Fisher matrix. In Sec. VI B we explore 



the possible parameter biases due to the uncertain cali- 
bration of the nonlinear power spectrum. 




k (h Mpc ) 



FIG. 4: The matter power spectrum at z — for three 
pairs of (Qx,w). Linear power spectrum correspond- 
ing to the fiducial spectrum is shown by the thin solid 
curve. Vertical lines delimit the interval which con- 
tributes significantly to the WL convergence power spec- 
trum, roughly corresponding to 100 < I < 10000. It 
can be seen that the ability to determine cosmological 
parameters will depend critically upon the knowledge of 
the nonlinear power spectrum. 



Figure ^ shows the matter power spectrum at z = 
for three pairs of (Qm,w). When fix or w are var- 
ied, the growth and normalization change affecting all 
scales equally. Varying fix also changes the transfer 
function (at k > 0.02 h _1 Mpc). On smaller scales 
(k > 0.2 h _1 Mpc), the non-linear power spectrum is fur- 
ther affected by dark energy. 



^ + l)P r /( 27r) oc|(l--^ 



ki 



(27) 



for I < fcir(z s ), and zero for I > kir(z s ). The plot of 
the convergence power spectrum is given in Fig. [5] for 
two values of k%. The multipole power peaks at / = 
3/5 kir(z s ). Assuming a survey with z s = 1, the scale 
at which the non-linear effects become significant, k rs 
0.2 h -1 Mpc, corresponds to I w 300. Our constraints 
mostly come from angular scales I ~ 1000, corresponding 
to k ~ 1 h _1 Mpc. The bulk of WL constraints therefore 
comes from non-linear scales. 




10000 



FIG. 5: Power spectrum of the convergence assuming 
matter power spectrum is a delta- function at k\, shown 
for two different values of k% . This shows the correspon- 
dence between physical and angular scales (for z s = 1 
and our fiducial ACDM cosmology). 



CONSTRAINTS ON DARK ENERGY 



C Angular Scale — Physical Scale Correspondence 

To illustrate the correspondence between wavenumbers 
k and multipoles I, let us assume the matter power spec- 
trum peaked at a single multipole k\ 

A 2 (k) = {S 2 )ki6(k-k 1 ) (26) 

normalized so that J A 2 (k)d\nk = (<5 2 ) (here (5 2 ) is 
the auto-correlation function of density contrast in real 
space) . Then, assuming for simplicity that all sources are 
at a single redshift z s , we have 



A The Fisher matrix formalism 

The fact that the relatively featureless depends 
upon a number of cosmological parameters directly leads 
to parameter degeneracies and limits the power of weak 
lensing to measure these parameters independently of 
other probes, even for the case of a full-sky survey. To 
estimate how accurately cosmological parameters can be 
measured, we use the Fisher matrix formalism E3. This 
method has already been used to forecast the expected 
accuracies from CMB surveys @, @, SNe la §, | 
and number counts |48[ and was found to agree very well 
with direct Monte-Carlo error estimation. Its consider- 
able advantage over Monte Carlo is that it does not re- 
quire simulations and analyses of data sets, but only a 
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single evaluation of a simple analytic expression. Fur- 
thermore, the Fisher matrix formalism allows easy inclu- 
sion of Bayesian priors and constraints from other meth- 
ods. 

The Fisher Matrix is defined as the second derivative 
of the negative log-likelihood function 



Fij = - 



d 2 InL 



dpidp 



3 I x 



E 



1 dPf dPf 
(AP*) 2 d Pl dp, 



(28) 
(29) 



where L is the likelihood of the observed data set x given 
the parameters p% . . . p n . The second line follows by as- 
suming that L is Gaussian in the observable Pf, which 
is a good assumption for small departures around the 
maximum. 

In practice we do not estimate the power spectrum at 
every multipole I, but rather bin Pf in 17 bins. We ex- 
plicitly checked that binning makes no significant differ- 
ence in our results (this is not surprising, as the conver- 
gence power spectrum doesn't have features that would 
get washed out by moderate-resolution binning). We con- 
sidered Pf at 100 < I < 10000, corresponding to angles 
between 1 arcminute and 2 degrees on the sky. Variations 
in the minimum and maximum I do not change any of our 
results, as very large and very small scales are dominated 
by cosmic variance and Poisson noise respectively. 

Finally, we need to choose steps in parameter direc- 
tions when taking numerical derivatives. We choose the 
steps to be 5% of the parameter values, making sure to 
take two-sided derivatives. 



B The fiducial cosmology and fiducial survey 

Finally, we need to choose the fiducial survey, i.e. sky 
coverage and depth of the survey. We do not consider any 
single experiment in particular, but rather adopt num- 
bers roughly consistent with proposed dedicated wide- 
field surveys expected to become operational in several 
years. We assume a survey covering 1000 sq. deg. down 
to a limiting magnitude R — 27; dependence of th e rc- 
sults upon these two parameters is discussed in Sec. VF. 



Surveys of this power are not yet operational, but are 
expected in the near future with results perhaps by the 
end of this decade. To convert from magnitudes to sur- 
face density of galaxies, we use the correspondence from 
Herschel and Hubble Deep Fields [Q , which for our fidu- 
cial numbers implies 165 gal/arcmin 2 . We assume that 
the only sources of noise are statistical: cosmic variance 
which dominates on large scales, and shot-noise domi- 
nant on small scales. We discuss the effect of systematics 
in Sec. |VI. 



C Parameter space 



Power spectrum of the convergence depends on 7 pa- 
rameters: fix, w, flhih , fish 2 , Sh, ti, and m„, where 
fls is the energy density in baryons (relative to criti- 
cal), n is the spectral index of scalar perturbations, and 
m v the neutrino mass summed over all species. In addi- 
tion, Pf depends upon the redshift distribution of source 
galaxies. Throughout, we use a fiducial model that fits 
well all experiments so far: fix = 1 — fin = 0.7 (flat 
universe assumed), h = 0.65, fish 2 = 0.019, n = 1.0, 
and Sh infer red from COBE measurements as described 
in Sec. IV B. The mass of neutrino species is quite uncer- 
tain, but, according to solar neutrino experiments, likely 
to be between zero and a few eV; we adopt m v = O.leV. 

We would like to get an insight in parameter degenera- 
cies, in particular between the equation of state ratio w 
and other parameters. To do that, we compute the corre- 
lation between w and other parameters. The correlation 
coefficient is given by 



p{w,Pi) 



Cov (w,Pi) 



y/Cov(w, w)Cov(pi,Pi) 



(30) 



where Cov(pi,pj) — is an element of the covariance 
matrix. Because imposing priors would alter the covari- 
ance matrix and confuse its interpretation, at this point 
we add no priors except for COBE normalization (10% 
in 5h) and perfect knowledge of galaxy redshifts. 

The most significant correlations are p(w, fix) — 
-0.96, p(w,fl B h 2 ) = -0.83, and p{w,m v ) = 0.81. We 
find that these and other correlations are typically very 
dependent on the fiducial model and the assumed prior. 
Finally, we examine the eigenvalues and eigenvectors of 
the Fisher matrix. The combination 0.61 fix + 0.21 w 
is determined to an accuracy of about 0.03; this is the 
best-determined combination containing significant com- 
ponents in fix and w directions. The most poorly com- 
bination of all is one almost entirely in the w-direction: 
0.97 w — 0.21 fix; it is determined to about 0.4. 



D Bayesian Priors 

Without any prior information on cosmological param- 
eters, weak lensing imposes very weak constraints on dark 
energy (and other parameters as well). The reason is that 
the power spectrum of the convergence is featureless, ow- 
ing to the fact that it represents the radial projection 
of the density contrast. Unlike the CMB spectrum, it 
lacks bumps and wiggles that would help break parame- 
ter degeneracies. Constraints rapidly improve, however, 
if the redshift distribution of source galaxies is known. 
We assume this to be the case; indeed, photometric red- 
shift techniques already show that distribution of source 
galaxies in weak lensing surveys will be determined inde- 
pendently of cosmological parameters (e.g., ]5Q] ) . Exact 
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knowledge of the source distribution is obviousl y a strong 
and perhaps unrealistic assumption, and in Sec. [Vj we ex- 
plore what happens when the uncertainties are included. 

There is no reason to expect that any cosmological 
probe alone should carry the burden of determining all 
cosmological parameters. Indeed, a number of cosmolog- 
ical parameters are already pinned down quite accurately 
by other means. In about 10 years, when powerful weak 
lensing surveys we consider complete their observational 
programs, parameters inferred from the CMB (such as 
Q,Mh 2 , fish 2 and n) will be determined to an accuracy 
of several percent [Ba. The neutrino mass, on the other 
hand, is poorly known today, but in the near future it is 
likely to be constrained by a combination of CMB, Ly- 
a forest ]5]]] , as well as solar and atmospheric neutrino 
measurements. 

For these reasons, we include Gaussian priors on cos- 
mological parameters (other than tlx and w). We con- 
sider two sets of priors, and call them "COBE+photo-z" 
and "Planck (T)" . The former set of priors is a weak one: 
we only include the 10% uncertainty in COBE normal- 
ization and, as mentioned above, knowledge of the distri- 
bution of background galaxies. The latter set is a mod- 
erate one, corresponding to the COBE+photo-z prior, 
plus the constraints expected from the Planck mission 
with temperature information only (Table 2 of Ref. [[46||): 
a(\ntl M h 2 ) = 0.064, a(lntl B h 2 ) = 0.035, a(n) = 0.04, 
and tj(m v ) = 0.58[f] . We note, however, that details 
of the second prior do not change the results much; for 
example, using the considerably weaker assumptions cor- 
responding to MAP mission (with temperature only) in- 
stead of Planck (T), errors in tlx and w degrade by only 
10% and 5% respectively. Similarly, using the very strong 
prior of Planck constraints (temperature and polariza- 
tion) combined with those from Sloan Digital Sky Survey 
(SDSS[|] ), the constraints improve only by about 20%. 
The reason for this weak dependence on the prior is easy 
to understand: by assuming the knowledge of the dis- 
tribution of source galaxies and adding other priors, we 
have broken the major degeneracy between tlx, w and 
other parameters; further information on other param- 
eters leads to small improvements in the constraints on 
dark energy. 



[f] Strictly speaking, the correct way to add the CMB priors would 
be to add the WL and CMB Fisher matrices. This procedure would 
correctly account for breaking of the WL parameter degeneracies by 
the CMB. We opt, however, to just add the priors to the diagonal 
elements of the WL Fisher matrix. This effectively assumes other 
parameters to be constrained within some limits, regardless of what 
experiment those constraints come from. 

[|] http: / / www.sdss.org 
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FIG. 6: 68% C.L. constraints on VLm and w for two val- 
ues of m v . We assume a survey of 1000 sq. deg. down to 
27th magnitude in i?-band, assume knowledge of the dis- 
tribution of source galaxies, and impose the Planck (T) 
prior on other parameters [Em . The strength of the con- 
straints does not depend sensitively on the set of priors, 
but does depend on the fiducial model (e.g., the neutrino 
mass). For orientation, current 1-tr constraints from 42 
type la supernovae 0] are also shown. 



E Results 

An example of the constraints that weak lensing can 
impose on dark energy is shown in Fig. ^| Here we show 
the 68% constraint regions for our fiducial WL survey 
(1000 sq. deg. down to 27th mag) with several sets of 
priors on other parameters. The ellipse is oriented so 
that increase in w is degenerate with increase in Hm, 
which is opposite of what we would expect; this is due to 
the fact that we assume galaxy redshifts to be known [§] 
. Table 1 lists the uncertainties using two sets of priors. 
Weak lensing is potentially a strong probe of dark energy: 
the 1-ct uncertainties in tlx and w are ~ 5% and 20-40% 
respectively (depending on the set of priors), which is 
somewhat weaker than statistical errors expected from 
future SNe la and number-count surveys. We emphasize 
that these numbers are the best ones possible given the 
survey specifications; systematic errors may weaken the 
constraints (see Sec. VI). It is also true, however, that 
weak lensing tomography can sig nifica ntly improve these 
constraints (more on that in Sec. VF ). 



[§]ln general, priors on other cosmological parameters will change 
the orientation of the constraints in the Qm~ w plane. 
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TABLE 1 
Constraints on dark energy 

Prior 

COBE + photo-z Planck (T) 



a(w) 



0.08 
0.36 



0.04 
0.19 



Top panel of Fig. ?] shows the dependence of the un- 
certainties in fix and w on the sky coverage, holding the 
depth of the survey fixed at 27th mag. Lower panel of 
the same Figure shows dependence of the uncertainties on 
the depth of the survey, holding the sky coverage fixed at 
1000 sq. deg. The constraints on w depend quite strongly 
on the depth of the survey — for example, constraint on 
w would improve by a factor of two by increasing the 
coverage of the survey to 5000 sq. deg. The dependence 
on the depth is also significant, but probably complicated 
by some practical problems; for example, galaxy overlap. 
Therefore, future surveys with very deep and/or wide sky 
coverage will be especially effective probes of dark energy. 



F Power spectrum Tomography 

One way to extract more information out of the data 
would be to divide the lensed galaxies in several red- 
shift bins and measure the convergence power spectrum 
in each bin, as well as the cross power spectrum between 
bins. This procedure, the power spectrum tomography, 
should be fully feasible with upcoming surveys because 
redshifts of source galaxies are going to be known quite 
accurately through photometric techniques. 

Following the formalism of Hu we compute the 
parameter constraints when source galaxies are separated 
in redshift. Of the several slicings in two bins we tried, 
the most effective division was below and above z = 1.0 
(Fig. ||, top panel). In this case, the constraints on fix 
and w improve by a factor of 3 and 1.4 respectively, for a 
Planck(T) prior (Fig. |). For the weaker COBE+photo-z 
prior, the improvement is even more significant: a factor 
of 5 and 3 improvements on fix and w respectively. 

We also consider an optimistic scenario where galax- 
ies can be separated in 10 redshift bins (Fig. [|, bottom 
panel). Whether or not and how accurately something 
like this can be done using photometric redshift tech- 
niques is presently under investigation (Eiscnstein, Hu 
and Huterer, in preparation). The constraints on Qm 
(or fix) and w further improve: o~(flx) = 0.012 and 
a(w) = 0.07 (Fig. §). 

Subdividing the galaxy population in more than two 
redshift bins leads to fairly limited improvements in pa- 
rameter determination; this is due to high correlations 
(~80%) between the power spectra in different bins [ p2[ . 
Nevertheless, tomography clearly adds valuable informa- 
tion on cosmological parameters and should be pursued 
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FIG. 7: Dependence of a(flx) and a(w) on the survey 
parameters. In each case we assume Planck (T) priors 
on other cosmological parameters. Diamonds denote the 
fiducial values. Top panel: 1-a uncertainties on fix and w 
as a function of sky coverage of the survey. We assume a 
fixed depth of 27th magnitude in i?-band. Bottom panel: 
1-a uncertainties on fix and w as a function of depth of 
the survey, assuming a fixed sky coverage of 1000 sq. deg. 



with data from future WL experiments. In order to accu- 
rately assess and optimize this technique, further study 
considering realistic accuracy of photometric redshifts is 
necessary. Using simplified assumptions (in particular, 
no "leakage" of galaxies between bins), we have shown 
here that separation of galaxies in redshift easily leads to 
a factor of a few improvement in measuring fix and w. 

We now discuss whether a specific signature of certain 
dark-energy models can be detected with WL surveys. 



G Detecting the Dark-Energy Clustering? 

Evolving scalar fields, or quintessence, are a particular 
class of candidates for dark energy (e.g. |33], [sBf). 
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FIG. 8: The divisions of source galaxies in rcdshift we 
used in order to implement the tomography. Top panel: 
A simple division in two redshift bins. Bottom panel: A 
division in 10 redshift bins. 
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FIG. 9: The improvement in the constraint on U,m and 
w due to tomography. The 68% C.L. constraint regions 
correspond to 1, 2 and 10 divisions in redshift (largest to 
smallest ellipse), and are all computed using the Planck 
(T) prior. 



One signature of quintessence is that it generally clusters 
around and above the Hubble radius scale Hq 1 . We ask: 
is it possible to detect this clustering in wide-field weak 
lensing surveys? 

The clustering of quintessence is reflected in the in- 
crease in the transfer function on very large scales. The 
effect is more pronounced for larger fix and larger w, 
and explicit forms for Tq(]c, z) are given in Refs. J35[ and 
|56f ; here Tq is the transfer function that takes clustering 
into account. 




FIG. 10: The effect of clustering of quintessence on the 
convergence power spectrum for a fiducial cquation-of- 
state ratio w — —1/3. The error bars correspond to 
the cosmic variance for a full-sky weak lensing survey. 
Clustering affects the I = 1 multipole the most, but even 
there the effect is buried within cosmic variance. 



Clustering changes the matter power spectrum on large 
scales, which in turn alters the convergence power spec- 
trum at lowest multipoles. In Figure |l0| we show an 
optimistic scenario[f] with w — —1/3 with and with- 
out clustering taken into account. We used exact formu- 
lae for the convergence power spectrum [ pT| , since Lim- 
ber's approximation breaks down at lowest multipoles. 
Even though the effect on the matter power spectrum 
is significant (T<g(fc, z) /T\(k, z) 2.0 at k ~ Ho in this 
case), the convergence power spectrum changes notice- 
ably only at I = 1, and even there only by ~ 30%. As 
this Figure shows, the effect is buried deeply in the cos- 
mic variance even for a full-sky WL survey. Therefore, 
it is unlikely that WL alone can detect the clustering 
of quintessence. However, cross-correlation of WL and 
other methods (e.g. the CMB) may be more promising; 
see Refs. § ||. 



[t]In a sense that a more positive w leads to more clustering. 
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VI. SYSTEM ATICS AND BIASES 

A Observational issues 

It is difficult to overemphasize the importance of con- 
trolling the various systematic errors that generically 
creep into the WL observing process. These include 
shear recovery issues, anisotropic point-spread function, 
the quality of seeing, and instrumental noise (for a nice 
study of systematic effects, see Ref. 60 ). There is also 
the effect of overlapping galaxies, which is expected to 
be especially pronounced in very deep surveys, but might 
be overcome using the photometric redshift information 
(M. Joffre, private communication). Finally, the ob- 
served galaxies might be intrinsically aligned due to cou- 
pling of their angular momenta or a similar mechanism 
([|Gl| p2| , |63| and references therein); this has already been 
observed |p4 |65| . These effects may masquerade as the 
signal itself, and make the extraction of ellipticity corre- 
lations very difficult. In our analysis, we have assumed 
that these problems will be resolved, and that the dom- 
inant uncertainty will be the cosmic variance on large 
scales and Poisson noise on small scales. In that sense, 
our results (for a given parameter space, set of priors, 
and fiducial survey strategy) may be optimistic. On the 
other hand, rapid advances in our understanding of weak 
lcnsing techniques, as well as the prospects of powerful 
future surveys, indicate that in a few years we can ex- 
pect a much better understanding of the aforementioned 
problems. 



B Dependence upon nonlinear power spectrum 
and galaxy distribution 

In addition to observational systematics that need 
to be controlled, theoretical prediction for the angular 
power spectrum of the convergence is also uncertain. 
Uncertainties in the nonlinear matter power spectrum 
(NLPS) and in the redshift distribution of galaxies are 
especially significant, as they are difficult to quantify and 
were not included in our analysis. We now discuss these 
two ingredients in more detail. 

As can be seen from Figs. || and ||, most of our con- 
straints come from nonlinear scales. Therefore, knowl- 
edge of the NLPS is crucial in order to compare exper- 
imental results with theory. However, this quantity is 
perhaps the most uncertain ingredient in the prediction 
for the power spectrum of the convergence. The NLPS is 
traditionally obtained by running N-body simulations for 
several cosmological models and deriving a fitting func- 
tion to the simulated nonlinear power spectra. For the 
models with dark energy we consider, either PD or Ma 
QCDM fitting functions can be used. The latter was cal- 
ibrated for quintessence models in a flat universe, and 
tested at w = —2/3, —1/2 and —1/3 and Qm — 0.4 and 
0.6. 



Even with this solution, the intrinsic uncertainty of 5- 
15% in the NLPS is significant (recall, the transfer and 
growth functions are accurate to just a few percent). 
To illustrate the importance of knowing the NLPS ac- 
curately, let us for the moment assume that the true 
NLPS at w = —1 is that given by the formula of PD. 
Let us further assume that, not knowing this, we adopt 
the Ma QCDM prescription to compute the theoretical 
power spectra. We now compute the bias in cosmological 
parameters due to this "erroneous" assumption. Let us 
write the cosmological parameter values as 



Pi = pi + 5p t 



(31) 



where pi is the true value, pi the measured value, and Spi 
the bias due to using the "wrong" NLPS. Assuming that 
these biases are small, it is easy to show that p6| 



(AP«)2 



r)P K 
dpj 



(32) 



where is the ubiquitous Fisher matrix, Pf (Pf) is the 
"erroneous" ("true") power spectrum, and sum over j is 
implied. The results of this exercise are given in Table 
2 where we consider our fiducial survey with Planck (T) 
prior. The biases in fix and w are 2.5 and 4.8 times the 
1-cr uncertainties in these parameters! Even though these 
numbers may not be accurate because the approximation 
6p <C p necessary to use Eq. (|||) obviously did not hold, 
one can still conclude that the biases are very significant. 
Therefore, we need a more accurate knowledge of the 
NLPS. 

TABLE 2 
Parameter Biases 



Due to 


"wrong" NLPS 


Due to 


"wrong" n(z) 


Pi |bias 


bias \/a(pi) 


|bias| 


bias \/a(pi) 


fix 0.09 


2.5 


0.04 


1.2 


w 0.92 


4.8 


0.57 


3.0 



Fortunately, the NLPS obstacle is surmountable. It is 
a matter of running powerful N-body simulations that 
include dark energy, on a fine grid in w (and ra„ and 
other parameters, if necessary). Because we are only in- 
terested in the matter power spectrum (not galaxy power 
spectrum, which includes bias), N-body simulations can 
in principle give the NLPS to a very high accuracy. Once 
this is achieved, weak lensing will regain much of its 
power to probe dark energy. 

Another quantity that may not be known to an ex- 
tremely high accuracy (although we assumed so) is the 
redshift distribution of source galaxies n(z). Indeed, cur- 
rent photometric redshift techniques can determine red- 
shifts to an accuracy of ~ 0.1, depending on the redshift 
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(e.g., which leaves room for error, both statistical 

and systematic. To include the uncertainty in n(z), some 
authors (e.g., [jl6[ |2^, parameterized the redshift dis- 
tribution by one parameter only. However, the realistic 
uncertainty in n(z) is much more difficult to quantify. To 
assess the effect of an uncertainty in the redshift distri- 
bution, we assume that the true distribution is given by 
Eq. (||) with zq = 0.5, while we "erroneously" assume the 
same form with zq — 0.55 (recall, n(z) peaks at z = 2zo). 
The biases in fix and w are given in Table 2, and are 1.2 
and 3.0 times the unbiased 1-a uncertainties in these pa- 
rameters, respectively. Just as in the case of the NLPS, 
we conclude that accurate knowledge of the redshift dis- 
tribution of galaxies will be crucial if weak lensing is to 
achieve its full potential. 



C Power spectrum covariance 

Yet another important issue that we ignored so far is 
covariance of the convergence power spectrum. The shear 
(or convergence) field is expected to be non-Gaussian due 
to nonlinear gravitational processes. Therefore, measure- 
ments of Pf are generally going to be correlated, imply- 
ing a non-zero four-point function (or its Fourier ana- 
logue, the trispectrum). The covariance will be especially 
pronounced at high multipoles. For a survey down to a 
limiting magnitude of R ~ 25, the effect of power spec- 
trum covariance appears to be small: Cooray & Hu Q 
have used the dark-matter halo approach to compute the 
power spectrum as well as the trispectrum, and found 
that the non-Gaussianity increases errors on cosmolog- 
ical parameters by about 15%. Although this effect is 
small enough to be ignored with current datasets, it will 
be important to take it into account when interpreting 
results from upcoming deep surveys because the covari- 
ance on small scales is likely to significantly degrade the 
cosmological constraints. Restricting our analysis (with 
COBE+photo-z prior) to multipoles / < 3000 degrades 
the constraints on fix and hi by a factor of 5. Clearly, 
information from small scales is important, and it will be 
necessary to carefully assess the impact of power spec- 
trum covariance for deep WL surveys. 



VII. THREE-POINT STATISTICS AND DARK 
ENERGY 

We now turn to three-point statistics of the weak lens- 
ing convergence. Unlike the CMB temperature fluctu- 
ations which may or may not be Gaussian, weak lens- 
ing convergence almost certainly does not obey Gaussian 
statistics. In this Section, we illustrate the dependence of 
the bispectrum and skewness of the convergence on dark 
energy, and show that they present a promising avenue 
that can lead to the dark component. 



A Preliminaries 

The bispectrum of the convergence B£j 2 j 3 is defined 
through the three-point correlation function of the con- 
vergence in multipole space 



h h h i g K 



mi m 2 TO3 
and can further be written as 



(33) 
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(2Zi + l)(2Z 2 + l)(2Z 3 + l) 
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r(X) 



1-2 



Ji 

r (x) ' r(x) ' r(x) 



■X 



(34) 



The bispectrum is defined only if the following rela- 
tions are satisfied: \lj — lk\ < h < \lj + h\ for {i, j, k} £ 
{1,2,3} and l\ + l 2 + ^3 is even. The term in paren- 
theses is the Wigner 3j symbol, which is closely related 
to Clebsch-Gordan coefficients from quantum mechan- 




|6SJ). W( X ) is the 
To compute the bis- 



ics (for its properties, see Refs 
weight function defined in Sec. 
pectrum of the convergence, therefore, we need to supply 
the matter bispectrum B(k\, k 2l k^,z). The latter quan- 
tity can be calculated in linear theory (that is, on large 
scales), but, just as in the case of the matter power spec- 
trum, it needs to be calibrated from N-body simulations 
on nonlinear scales. Here we adopt the fitting formulae 
of Scoccimarro & Couchman (|JFo|; heretofore SC) which 
are based on numerical simulations due to VIRGO col- 
laboration |7])| . The matter bispectrum is defined only 
for closed-triangle configurations (k\ + k 2 + k% = 0) and 
is given by 



B(k 1 ,k 2 ,k 3 )^2F 2 {k ll k 2 )P(k 1 )P(k 2 ) + cyc\. (35) 
where P(k) is the matter power spectrum and 



F 2 (ki,k 2 ) = -a(n, k x ) a(n, k 2 ) + 



2 k x -k 2 



7 \ fci k 2 



c(n, fci) c(n, k 2 ) 



(36) 



n = <iln P/ d\n k, and functions a, b and c are given in 
SC. Although not explicitly tested on models involving 
dark energy, the fitting formula depends on cosmology 
only through the matter power spectrum; this weak de- 
pendence on cosmology is also borne out in high-order 
perturbation theory 74 . Therefore, we decide to use the 
SC formula to illustrate the dependence of three-point 
statistics on dark energy. 
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FIG. 11: The quantity I 2 Bj/ 2 / (2ir) , involving 
equilateral-triangle configurations of the bispectrum in 
multipole space. We use this quantity to illustrate how 
the bispectrum depends on dark energy. The variance in 
Bui is roughly two orders of magnitude larger than the 
signal. 



In Fig. [ll] we show the quantity I 2 y/Bfc / '(2tt) fr| for 
w = —1 and w = —0.5; here B^ n is the equilateral tri- 
angle configuration of the bispectrum [f] . Since roughly 
B oc P 2 and B has little other dependence on dark en- 
ergy, we expect that I 2 Bf a / [2n) varies with w simi- 
larly as P — and this is correct (compare Figs. || and 
|ll|) . Therefore, the bispectrum appears to be an ex- 
cellent probe of dark energy. Things are complicated, 
however, by the large cosmic variance of a bispectrum. 
Although computing variance of B involves a daunting 
task of evaluating the six-point correlation function of the 
convergence, this quantity can be computed under an as- 
sumption of small departures from Gaussianity |75|, [76[ . 
For the equilateral triangle configuration of the bispec- 
trum we show, this estimate indicates that the cosmic 
variance is about two orders of magnitude larger than 
the bispectrum signal itself, roughly independently of I. 
Therefore, it is unlikely that a single configuration of the 
bispectrum can be used to probe dark energy. However, 
one should be able to find an optimal combination of 
configurations in order to maximize the amount of infor- 
mation. We relegate this problem to future work. 

Next we discuss the dependence of skewness on dark 
energy. Skewness is defined as 



Ss(0) 



(n 2 (9)) 2 



(37) 



where 



( K *tf)) = J_ J- ■ /( 2 ^ + 1 )( 2/ 2 + 1 '(2/, I.) 



(38) 



hhh 



4tt 



x(o o '5)w.( s )^(W.( J )- 

(39) 

are the second and third moments of the map smoothed 
over some angle theta, and W; (9) is the Fourier transform 
of the top-hat function: Wi{9) = 2J 1 (W)/(W). Skewness 
effectively combines many different bispectrum configu- 
rations, and its variance should be much smaller than 
that of Bi ± i 2 i 3 . Its disadvantage is that measurements on 
different scales are correlated. 
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FIG. 12: Skewness of the convergence for two values of 
w. Error bars are from simulations by White & Hu |77j| 
on scales they explore and for a field of 36 sq. deg. 



Fig. |l2] shows skewness for two values of w. Roughly 
speaking, S 3 cx B K /P 2 , and although B K and P 2 both 
decrease with increasing w, the P 2 term prevails — hence 
the scaling of 53 with w. The error bars shown are those 
from White & Hu []77j for their WL simulations corre- 
sponding to the ACDM model, and for a field of 36 sq. 
deg. Although the dependence of skewness on dark en- 
ergy is significant, there are several obstacles. As in the 
case of the matter power spectrum, the fitting formula 
for the bispectrum is accurate only to about 15% (rms 
deviation) for ACDM models and not yet calibrated for 
dark energy models. More seriously, the measurements 
of skewness are likely to be highly correlated — in fact, 
van Waerbeke et al. (72| find that correlation between 
skewness measurements (for the top-hat filter we use) is 
close to 100%. 



[f]We set m„ = in this Section. 
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In conclusion, our preliminary analysis indicates that 
the three-point statistics of the weak lensing conver- 
gence are sensitive to the presence of dark energy, mainly 
through the dependence of the matter power spectrum. 
More work is needed, however, in order for the three- 
point statistic to become an effective probe of the missing 
component. This will include sharpening the predictions 
for the three-point function in the nonlinear regime, and 
finding optimal configurations of the bispectrum to probe 
dark energy. 

VIII. DISCUSSION AND CONCLUSIONS 

Recent results coming from type la supernovae, CMB, 
and large-scale structure surveys make a strong case for 
the existence of dark energy. It is therefore important to 
explore how upcoming and future surveys can be used to 
probe this component. In this work, we explore the power 
of weak gravitational lensing to probe dark energy via its 
measurements of the power spectrum of the convergence. 

Dark energy modifies the convergence power spectrum 
by altering the distance-redshift relation, as well as the 
matter power spectrum. The dependence on dark en- 
ergy is therefore somewhat indirect, and cannot be easily 
disentangled from the effect of other parameters (£l]\ih 2 , 
fish 2 , n, m v ). Because of this, one would not expect WL 
to be an efficient way to probe dark energy. Neverthe- 
less, we find that with the proposed very wide and deep 
surveys, WL can be an important probe of dark energy, 
on a par with SNe la and number counts. We consider 
a generic future survey covering 1000 sq. deg. down to 
a limited magnitude of R = 27, cosmic-variance limited 
on large scales and Poisson-noise limited on small scales. 
With photometric redshift information and constraints 
on other parameters that would be expected from the 
Planck experiment with temperature information only, 
we find that such a survey is able to constrain Qx and 
w to between a few percent and a few tens of percent, 
depending on the fiducial model and a chosen set of pri- 
ors. The constraints are in general stronger for wider and 
deeper surveys, and depend on the fiducial model (e.g., 
the neutrino mass). Accurate knowledge of the redshift 
distribution of source galaxies will be crucial; we find that 
an error of only 0.05 in the peak of the redshift distribu- 
tion can bias the results. 

There are important caveats to this result, however. 
Most information from WL comes from nonlinear scales, 
where the evolution of density perturbations is difficult 
to track analytically and understood mostly through N- 
body simulations (restricting the analysis of P ; K only to 
linear scales with I < 100 would lead to extremely weak 
constraints on cosmological parameters due to cosmic 
variance). The nonlinearities potentially lead to at least 
two sources of systematic error. First, the power spec- 
trum measurements are likely to be strongly corre- 
lated at multipoles of several thousand and higher. This 
is especially true for planned deep surveys (down to a 



limited magnitude of R ~ 27 or higher), and these cor- 
relations will likely degrade the constraints on Qx and 
w. 

Second, although the nonlinear power spectrum has 
been calibrated quite accurately for ACDM models, most 
notably through the PD formula, it remains poorly ex- 
plored for models with general equation of state w, mas- 
sive neutrinos, and significant baryon density. We explic- 
itly showed that the lack of knowledge of the dependence 
of the nonlinear power spectrum on w can easily bias 
the constraints on fix and w. Therefore, a better under- 
standing and calibration of the NLPS is absolutely crucial 
in order to use WL as a tool of precision cosmology. 

This problem can be turned around, however. One 
could use the weak lensing measurements (combined, per- 
haps, with information from CMB measurements, SNe 
la, and other probes) in order to constrain the nonlin- 
ear power spectrum. This constraint could be very in- 
teresting, given the strong dependence of the NLPS on 
cosmological parameters. 

Predictions for the three-point statistics of WL are 
quite uncertain at present, especially for models involving 
dark energy. This does not mean they will not become 
effective probes of the missing component in the future. 
We estimate the equilateral bispectrum configuration, as 
well as skewness, for two values of w and show that de- 
pendence on w is significant. Although these two quan- 
tities are plagued by large cosmic variance and highly 
correlated noise respectively, by clever choice of bispec- 
trum configurations one might be able to increase the 
signal-to-noise ratio and extract useful information on 
dark energy. 

There are other ways to use weak lensing as a probe of 
cosmology which we did not discuss. For example, one 
could use WL to identify clusters of galaxies at redshifts 
< z < 3 (M. Joffre et al., in preparation). Comparing 
the measured number density of clusters to the predic- 
tion given by the formalism of Press & Schechter |78[ ] 
gives constraints on cosmology. Another idea is to mea- 
sure the angular power spectrum of clusters (detected 
through WL) at different redshifts f7{|; this gives a direct 
measure of the angular diameter distance as a function 
of redshift. The advantage of this approach is that only 
the linear matter power spectrum is required; further- 
more, the mass function and profiles of clusters need not 
be known. These two methods will provide constraints 
complementary to those from the galaxy shear. 

Weak gravitational lensing is likely to provide a wealth 
of information not only on the matter distribution in the 
universe, but also on the amount and nature of dark en- 
ergy. We have considered the basic program of measur- 
ing the convergence power spectrum, and found that very 
wide and deep surveys could provide information comple- 
mentary and comparable to that from other cosmological 
probes. Other statistics (various bispectrum configura- 
tions, cross-correlation of WL and the CMB, etc.) are 
likely to further increase the power of weak lensing and 
make it an important probe of dark energy. 
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